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Abstract 

A four-nucleon calculation of the total ^He photodisintegration cross section 
is performed. The full final-state interaction is taken into account for the 
hrst time. This is achieved via the method of the Lorentz integral transform. 
Semi-realistic NN interactions are employed. Different from the known partial 
two-body ^He( 7 , n)^He and ^He( 7 ,p)^H cross sections our total cross section 
exhibits a pronounced giant resonance. Thus, in contrast to older ( 7 , np) 
data, we predict quite a strong contribution of the ( 7 , np) channel at the 
giant resonance peak energy. 
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The photodisintegration of “^He has received much attention in the last 25 years. Experi¬ 
mental work concentrated mainly on the two dominant two-body breakup channels (^He-fn, 
In a hrst round of experiments a rather strong peak of the giant dipole resonance 
was found, while more recent experiments hnd a much less pronounced peak. The suppres¬ 
sion of the two-body breakup peak was conhrmed in four-nucleon calculations that take 
into account the important hnal state interaction (ESI) via a semi-realistic NN potential 
imj^. Much less is known about the total a photoabsorption cross section (^He-fy X). 
In the vicinity and beyond the peak there are neither theoretical calculations that take into 
account ESI nor experimental total cross section measurements. 

The situation for the ^He photodisintegration seems to be settled only for the two-body 
breakup channels at lower energies. Yet the results are rather puzzling because it is not 
understood why the a-particle should have such a suppressed giant dipole resonance. Cross 
sections for transitions to other channels {''y,pn)d, ( 7 , 2p2n) and ( 7 , d)d obtained in the older 
experiment are very small and cannot influence the general picture at all. Furthermore, 
the new ( 7 ,p)^II and ( 7 ,n)^He data combined with those cross sections would lead to a 
bremsstrahlung weighted sum over the photoabsorption spectrum which is substantially 
lower than the well known model-independent sum rule estimate. A theoretical calculation 
of the total photoabsorption cross section would certainly help to get a better understanding 
of these problems, since the giant resonance is in principle a feature of the total cross section. 

In the present work the theoretical calculation of the total cross section is carried out 
with consideration of the full ESI. Previously the ESI was taken into account completely 
only below the three-body p + n + d breakup threshold E-y = 26.1 MeV [|,||. For the 
two-body breakup the resonating group calculation of Ref. [Q was extended to somewhat 
higher energies taking into account ESI due to other channels approximately. At > 50 
MeV the two-body reactions were treated in the plane-wave approximation [^. 

We calculate the total photoabsorption cross section in the whole energy range below 
the pion threshold. We consider the El transition in the long-wavelength limit using the 
unretarded dipole transition operator 
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i=l 

In this way we take into account meson exchange currents via the Siegert theorem. The E2 
contributions to the total cross section are small even at high photon energy [Q and they 
tend to cancel with the El retardation contributions |p. Our nuclear hamiltonian includes 
central even local NN potentials and the Coulomb interaction. 

We can write down the total photoabsorption cross section as 

atot{E^) = 4:n^{e^/hc)E^R{E^) , 

where R is the dipole response function, 

R{E,) = J df\{'^f\D,\'^o)\"S{Ef -Eo- E ,). 

Here Tq is the a-particle wave function and Tj are hnal state wave functions normalized as 
(T/|T/-/) = S{f — /'). In the above relations we neglect the very small nuclear recoil energy. 
We calculate the response function R via evaluation and subsequent inversion of its Lorentz 
integral transform, a method we proposed for the response of an arbitrary N particle system 
to an external probe 0. The method has already been successfully applied for obtaining 
the accurate longitudinal (e, e') response functions of the two-, three-, and four-nucleon 
systems [|]-|^. The transform C{(j) of the response R is found as 

/:(a) = (T(a)|T(a)), (1) 

T being the solution to the Schrodinger-like equation 

{H-Eo + a)^{a)=Q ( 2 ) 

with the source-term Q = The function T is localized and continuum calculations 

are thus avoided in our approach. 

We use the same NN potential model, Trento (TN) potential, as in our work on the 
longitudinal response function |P]. We also consider the Malfliet-Tjon (MT) I+III potential 
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which was used in Ref. p| for calculating the reaction 7 +'^He —^-^He +n. We use the 
value of A = 1.555 fm“^ entering the attractive part of the MT potential as listed in Ref. 

I^] . This value just leads to correct low-energy parameters of NN scattering as given in 
Ref. 0. In some "^He bound-state calculations the value A = 1.55 fm“^ listed in Ref. P has 
been used that leads to an increase in the i?6(^He) value by about 1.4 MeV. 

Our a-particle wave function Tq is an eigensolution for the same NN potential. The 
corresponding matter r.m.s. radia and binding energies are 1.41 fm and 30.5 MeV for the 
TN potential, and 1.43 fm and 29.2 MeV for the MT potential. The latter value is close to 
those reported in the literature, see Ref. p|. The binding energies are reasonable as compared 
to the experimental value of 28.3 MeV, and the radia are close to the experimental value of 
1.45 fm. 

In Fig. 1 we show the s-wave phase shifts (no Coulomb interaction included for ^S'o) of 
both potential models in comparison to those of a realistic interaction (Paris potential [pT[| ). 
It is evident that MT and TN potentials do not lead to signihcantly different phase shifts 
than the Paris potential. The scattering length equals to -17.9 fm for the TN potential 
and -23.3 fm for the MT potential, so the MT potential is a little bit more attractive in the 
^S'o channel than the TN potential. The TN scattering length is close to the value of nn and 
pp (no Coulomb force) scattering (ann(^5'o) = —17.6 fm for Paris potential), while the MT 
scattering length is close to that of np scattering (a„p(^S'o) = —23.7 fm). 

We solve Eq. (2) for L = T = 1 and S' = 0 with the help of the correlated hyperspherical 
expansion and the hyperradial expansion of the same form as in Ref. |P|. The K^ax value 
equals to 7. The a value in Eqs. (1), (2) is of the form + iai with aj = const, and the 
values of cx/ = 20 MeV and 5 MeV have been employed. In Fig. 2 the convergence of the 
transform, Eq. (1), with respect to K^ax is shown for u/ = 20 MeV for the MT potential. 
While inverting the transform the true low energy behavior [E^ — {E^)miri\^^‘^ have been 
incorporated into our trial response. The inversion has been performed both for aj = 20 
MeV and for a combination of the transforms with aj = 5 and 20 MeV chosen so that the 
former transform gives a predominant contribution to the very steeply rising low energy 
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wing of the response and the latter ones to its high energy wing. The responses obtained 
in these two versions practically coincide with each other. The transforms in Fig. 2 with 
Kmax = 5 and 7 lead to practically identical responses and that for Kmax = 3 is also not 
very different. For the TN potential one hnds a similarly good convergence in Kmax as well. 

Besides the checks of the convergence, the overall test of the hnal results is pro¬ 
vided by sum rule calculations. We compare the bremsstrahlung weighted sum at = 
J^h atot{E^)E~^dEj and the TRK sum aTRK = <Jtot{E^)dE^ = 59.74(1 -|- k) MeV mb 
calculated with our cross sections with an independent calculation of these quantities us¬ 
ing the sum rules ( /?ic){^o\D^D^\^o), n = (ToliHz, [V, T)^]]|\I/o)(m/h^)R/77Z). 

The sum rule values are ab = 2.41 mb, k = 0.727 for the TN potential and at = 2.48 
mb, K = 0.684 for the MT potential. By integrating our cross sections explicitly we obtain 
cTfe = 2.40 mb, K = 0.754 for the TN potential and at = 2.48 mb, k = 0.712 for the MT 
potential. The agreement of the at values with the sum rules is perfect that reflects a good 
accuracy of the low energy wings of the responses obtained. The resulting relative deviations 
from the TRK sum rule are about 1.5% for both potentials. 

One may note that the k, values for the potentials we use are lower than those provided 
by fully realistic NN interactions. The latter values range from 1.0 to 1.3 Ha, thus we 
underestimate aTRK by 15-25%. We believe that the main part of the missing strength 
should lead to an increase of the cross section at higher energies, while our potential models 
should provide quite realistic results up to the pion threshold. In fact a rough estimate of cJb, 
which we performed for realistic NN interactions, is close to the at values for our potentials. 
In any case, an increase in the k value would only strengthen our conclusions about the 
strong ( 7 , np) cross section which we predict below. 

At this point we should mention that our calculation is performed consistently with 
our semi-realistic hamiltonians, i.e. applying the Siegert theorem we use the energy eigen¬ 
values of the hamiltonian. However for comparison with experiment we perform the shift 
atot{E^) —>• atot{E^ -|- AEb), AEb being the difference of the calculated and experimental 
binding energies. In this way we obtain the proper breakup threshold thus correcting for 
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some overbinding of our a-particle. 

Unfortunately there are no direct experimental data on the ^He total photoabsorption 
cross section. Nevertheless we would like to make a comparison with experimental data. 
Therefore we proceed as follows. For the low-energy region we make interpolations of the 
( 7 ,n) data from and the ( 7 ,p) data from and sum up the resulting ( 7 ,p) and ( 7 ,n) 
cross sections (dotted curves in Figs. 3,4). Since the ( 7 ,d)d cross section can be safely 
neglected (see e.g. |^) this should lead to a rather good estimate for the total cross section 
below the three-body breakup threshold. Furthermore, we also show the cross sections 
of other low-energy experiments Assuming that ( 7 ,p) and ( 7 ,n) cross sections 

are more or less equal we double the experimental cross sections in order to have further 
estimates for the two-body breakup. Beyond 26.1 MeV they represent lower experimental 
bounds for the total cross section. In Fig. 3 these estimates are shown together with the 
calculated cross sections for MT and TN potentials. There is a rather good agreement of 
our responses with the estimated experimental two-body cross section up to the three-body 
breakup threshold. The MT potential leads to a slightly higher low-energy cross section 
than the TN potential that may be related to the somewhat stronger attraction in the NN 
^S'o channel. For the MT potential we hnd a similar agreement with experimental data as 
was found in for the same potential model for the ( 7 ,n)^He channel. 

Beyond the three-body breakup threshold our cross sections reveal further increase. 
Since theoretical as well as experimental results for the ( 7 ,p) and ( 7 ,n) cross sections show 
a flattening beyond the three-body threshold, the further increase has to be attributed to 
( 7 , np) reactions. Thus the ( 7 , np) channel increases the peak of the giant dipole resonance 
considerably. As can be seen in Fig. 4 it leads to a rather pronounced resonance peak. Also 
in the high-energy sector we show lower experimental bounds for the total cross section. 
They consist of the sum of the ( 7 ,p)^H and ( 7 ,? 7 ,)^He cross sections from Ref. p] and the 
doubled ( 7 ,p)^H data from Ref. 0 . From the comparison of these estimates with our 
theoretical total cross sections one would expect quite an important contribution of the 
( 7 , np) channel in the whole energy range. 
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Finally, we summarize our work. For the first time the total cross section of the a- 
particle photodisintegration was calculated in the framework of four-nucleon dynamics with 
full FSI. The results show a very pronounced peak of the giant dipole resonance. Therefore 
it seems that a typical many-body feature emerges also from a genuine few-body calculation 
of the four-nucleon system. The peak is considerably higher than the sum of the cross 
sections of the two important two-body breakup channels ^He-f-n). Thus we predict 

quite a strong contribution of the ( 7 , np) channel already at rather low energies. More 
experimental work is needed to conhrm this prediction. At present some data on ( 7 , np) 
with high statistics are available only beyond 80 MeV while the energy range between 
three-body breakup threshold and 80 MeV remains to be explored. 

The authors thank H.M. Hofmann for helpful correspondence. 
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FIGURES 

FIG. 1. NN scattering phase shifts of the partial waves (a) and (b) for the following 
potentials: TN (dashed curves), MT (dotted curves), and Paris (full curves). 

FIG. 2. The Lorentz transform for the MT potential with various K^ax values. 

FIG. 3. Theoretical results for the total ^He photoabsorption cross section at low energy 
with MT (dashed curve) and TN potentials (full curve). Also shown is the estimate for the 
two-body breakup (dotted curve with typical size of the experimental error), which is based on 
the experimental results of Refs. [15,16] as well as doubled experimental cross sections for {'y,p) 
[19] (open circles) and for ( 7 , n) [18] (triangles) and [20] (full circles) (for further explanation see 
text). The three-body breakup threshold is marked by an arrow. 

FIG. 4. As Fig. 3, but for an extended energy range up to 140 MeV. Estimate for lower exper¬ 
imental bound (dotted curve) and additional lower bound estimates with data from [3] (diamonds), 
[19] (open circles), and [21] (squares) (for further explanation see text). 
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